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Abstract 

Quantum estimation of the operators of a system is investigated by analyzing its 
Liouville space of operators. In this way it is possible to easily derive some general 
characterization for the sets of observables (i.e. the possible quorums) that are 
measured for the quantum estimation. In particular we analyze the reconstruction 
of operators of spin systems. 



1 Introduction 

Two fundamental restrictions limit the possibility of devising a state recon- 
struction method. On one hand, the quantum complementarity principle does 
not allow to recover the quantum state from measurements on a single sys- 
tem, unless we have some prior information on it. On the other hand, the 
no cloning theorem ensures that it is not possible to make exact copies of a 
quantum system, without having prior knowledge of its state. Hence, the only 
possibility for devising a state reconstruction procedure is to provide a mea- 
suring strategy that employs numerous identical (although unknown) copies 
of the system, so that different measurements may be performed on each of 
the copies. 

The problem of state estimation resorts essentially to estimating arbitrary 
operators of a quantum system by using the result of measurements of a set 
of observables. If this set of observables is sufficient to give full knowledge 
of the system state, then we define it a quorum. Notice that, in general, a 
system may allow various, different quorums. Quantum tomography was born 
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[1] as a state reconstruction technique in the optical domain, and has recently 
been extended [2] to a vast class of systems. By extension, we now denote 
as "Quantum Tomography" all unbiased quantum state reconstruction proce- 
dures, i.e. those procedures which are affected only by statistical errors that 
can be made arbitrarily small by increasing the number of measurements. To- 
mography makes use of the results of the quorum measurements in order to 
reconstruct the expectation value of arbitrary operators (even not observables) 
acting on the system Hilbert space. 

The purpose of this work is to present in a formally familiar manner (employing 
the Dirac notation also on operator space) a constructive method to derive 
tomographic formulas for quantum systems, at least for finite dimensional 
Hilbert spaces. This is achieved by giving conditions to build quorums and to 
check whether a given set of operators is a quorum. In this way, we obtain 
an extension of the recently proposed group tomography [2], where similar 
conditions were derived for systems with an underlying group structure. 

In Sect. 2 we give the definitions and the conditions to identify a quorum of 
operators by analyzing the space of operators of a system as a linear vector 
space. We derive a constructive algorithmic procedure to obtain tomographic 
formulas in the case of finite quorums. In Sect. 3 we give some examples of 
applications of the presented method in the domain of spin systems, where 
various different quorums are available [2,3]. 



2 General estimation 



Consider the set of system operators, i.e. the Liouville space C{7i). If we 
initially restrict ourselves to Hilbert-Schmidt operators in C(7i), then this set 
is itself a Hilbert space of operators, with the scalar product 



{A\B) = Tr 



"i t B 



It is then possible to employ all the properties of linear vector algebra, and 
to use the Dirac notation, by using the following definitions for bra and ket 
vectors: 



6 — 16) 

Tr[.6t] — <6| . . (2) 

In this vision, quantum tomography consists of expressing the operator A we 
want to evaluate as an expansion on the observables of the quorum as 
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\A) = J dx \C(x))(B(x)\A) , 
x 



(3) 



where \A) is a generic operator in £(H), \C(x)) (with x G X) is the set 
of quorum observables (C(x) is a generally complex function of a selfadjoint 
operator, hence it is observable in this sense), and the set (B(x)\ is the dual 
of the quorum. In ordinary notation, Eq. (2) is the tomography identity, i.e. 



A = dx Tr 



C{x) (4) 



A' 



Notice that the extension of the theory to non-normalizable vectors in the 
operator Hilbert space is immediate: one only has to require the existence of 
the trace of Eq. (4). If, for example, A is a trace-class operator, then we do 
not need to require B{x) to be of Hilbert-Schmidt class, since it is sufficient to 
require B(x) bounded. Through Eq. (3), the tomographic reconstruction pro- 
cedure is immediately obtained. In fact, by measuring the observables \C(x)) 
of the quorum, we can(Q) express the mean value of any operator (A) in terms 
of the eigenvalues of \C(x)) as 



(A) = dx 5>(m,x) A£> Tr[&(x)A] , (5) 
x 

where p(m, x) is the probability of obtaining the eigenvalue \ff when mea- 
suring the quorum observable C(x). 

Since we want Eq. (3) to be valid for a generic operator \A) in C(TC) [or also 
in a subspace of C(TC)), then we must require that the \C(x)) constitute a 
spanning set for the operator (sub)space, with the set of (B(x)\ acting as its 
dual. A spanning set is a generalized basis for a vector space: it is a complete 
set of vectors but it is not, in general, composed of linearly independent (or 
normalized) vectors. Define dual (B(x)\ of the set \C(x)) as the set constructed 
so to have 

(B(x)\C{x')) = Tt[B\x)C(x')} = 5{x,x') \fx,x' G X, (6) 
where 8{x,x') is a reproducing kernel for (B(x)\, i.e. 



J dxS(x,x')(B(x)\ = (B(x')\ . (7) 
x 

2 Eq. (5) is obtained by taking the expectation value of both members of Eq. (3) 
and by calculating the expectation value trace using the eigenvectors of the quorum 
observables \C(x)). 
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Since \C(x)) is a complete set, 5(x, x') is a reproducing kernel also for this set, 
i.e. 



J dx5(x,x')\C(x)) = \C{x')) . (8) 
x 

From linear vector algebra we obtain the following four equivalent definitions 
of spanning set: 

A set of vectors \C(x)) (with dual (B(x)\) is a spanning set 

i) V|i) G | A) = J dx \C(x))(B(x)\A), i.e. the tomographic identity, 

x 

namely Eq. (3). 

ii) \C n ) is complete, i.e. (no nonzero element is orthogonal to \C(x)) Vx): 

(ilC^r)) = (£(x)|A) = Vx G * ^ |i) = . (9) 
hi) the following operatorial identity resolution applies, 

jdx\C(x)}(B(x)\ = l, (10) 



where 1 is the identity super-operator, namely the operator acting on 
operators such that i[A] — A \/A G C(H). 
iv) j dx {A\C(x)}{B(x)\A} =\\ A\\ 2 = Ti[A^A] V\A}eC(H). 



x 



In the usual notation, these equivalent definitions write as [4]: 

i) A = J dx Tt[B\x)A]C{x). 

x 

ii) Tr[&(x)A] = Tr[A^C(x)} = 0Vx G X A = 0. 

hi) Jdx i,\C{x)mk\B\x)\l) = 5 zl S jk , where {|n» is a basis for the system 
x 

Hilbert space ri. 

iv) J dx Tr[A*C(x)]Tr[&{x)A] = Tr[iU] Vi G C(H). 
x 

In order to obtain the dual set (B(x)\ starting from a given set \C(x)), one in 
general has to solve the operatorial equation (6) that defines the quorum. For 
finite quorums, this resorts to a matrix inversion. An alternative procedure is 
now proposed. It derives from the Gram-Schmidt orthogonalization method 
[5], which allows to derive a basis starting from a complete set of vectors. 
Namely, one obtains a basis \yk), given the complete set \Ck) (assume for 
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simplicity that all \Ck) are non-zero and that in {|Cfc)} there are no couples 
of proportional vectors), recursively defined as 



W) = -k {\C t )-T. k j ^\v,)(vi\C t )) 

where N =|| |C ) || and N k =|| \C k ) - Ej=d II- Notice that in the 

recursion (11) one must take care of eliminating all the vectors \C k ) which are 
a linear combination of the with j < k. 

Write the identity resolution for the basis obtained with procedure (11), i.e. 



1 = E \Vk)(Vk\ = 



k=0 

k-1 



( 0) (Vol + E w I \Ck) - E l%-X%|Cfc> ] (Vkl (12) 



By using repeatedly Eq. (11) (expressing \yj) of Eq. (12) in terms of the |C n )s) 
and by reorganizing the terms in the sums, we can find the dual set (B n \ as 



^ 01 A^o NqNi { N N 2 Nq^Nz J Km 

, , (yij _ (jHgbKjfaj / (yi|C 3 > (yi|C 2 )fo|C 3 ) \ , , 
^ 11 iVi iViJV 2 I iViJV 3 NiN 2 N 3 J [ 1 



(13) 



Eq. (12) guarantees that it is possible to write 

i = £|c„><fl B |, ( 14 ) 

n 

which is just the definition iii [i.e. Eq. (10)] of spanning set. 

Summarizing, we described a method for deriving tomographic formulas for 
arbitrary systems. One must start from a set of operators C n he would like to 
use as a quorum, and verify that such a set is complete, i.e. that no nonzero 
element of C(H) is orthogonal to all C n : 



(A\C n ) = Trfitcy = Wn \A) = . (15) 

If the set is finite, then one can employ the orthogonalization procedure out- 
lined previously to derive the dual set. If the set is infinite discrete or con- 
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tinuous, then one can only resort to finding appropriate solutions for Eq. (6). 
Once the dual is known, the tomographic identity (3) can be written explicitly. 
The reconstruction procedure, in terms of the probabilities of measurements of 
quorum observables, follows straightforwardly and yields Eq. (5), which allows 
to obtain arbitrary operator expectation values in terms of quorum outcome 
probabilities. Of course, one may think of similar procedures based on different 
orthogonalization algorithms. 

Since no hypotheses were made on the structure of the system Hilbert space, 
the theory presented in this section is valid for any quantum system. In the 
following section we will give some example applications. 



3 Example of application: Spin Tomography 



Here we show an application of the theory presented in the previous section 
by rederiving the spin tomography [2,3], where various different quorums may 
be employed. 

The simplest possible example is a spin s — | system. In this case we expect 
that the Pauli matrix and the identity constitute a quorum (since any 2x2 
matrix can be written on such a basis). Take the quorum Q = {a x , a y ,a z ,l}: 
it is immediate to verify that it is complete. Since the quorum operators are 
orthogonal, i.e. a a ■ a a / = lS aa i (a, a' = x,y,z), using the Gram-Schmidt 
procedure it is immediate to obtain the dual set as C = {^cr x , \o y , |<t 2 , |l}. 
The expansion (3) of a matrix A is, thus 



E |*a><*i|A> + |i><i|i) 



(16) 



la=x,y,z 

which immediately yields the reconstruction procedure 



<^>= E E p(rn,H a )m Tr 



m=-i a=x <v> z 



Aa n 



2 



A 



(17) 



where p(m,n a ) is the probability to obtain the eigenvalue m = ±| while 
measuring S • n a . This equation allows the reconstruction of the expectation 
value of any spin s — \ operator A from the measurement of the spin in the 
x, y, z directions. 



For an arbitrary spin s, a possible quorum is given by the spin component in 
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all directions, i.e. the observable S-n (S being the spin operator and n a vector 
on the unit sphere). In order to find the dual (B\, consider the exponential 
of the quorum, i.e. D(ip,n) = exp(iifjS ■ ft), which satisfies definition iii [i.e. 
Eq. (10)] of spanning set. In fact, D(tp,n) constitutes a unitary irreducible 
representation of the group SU(2). The orthogonality relation between the 
matrix elements of the group representation D{g) of dimension d writes as [6] 



J dg D jr (g)Dl k (g) = -S jk 8 tr , (18) 

R 

where dg is the group Haar invariant measure, and V — J R dg. For SU(2), with 
the 2s + 1 dimension unitary irreducible representation D(ip, n), Haar's invari- 
ant measure is sin 2 ^sint? d-d d(p dip, and V = Arc 2 . Thus, the orthogonality 
relations in this case are given by 



-^-t- Jdnjdij sin 2 |(j|e^- s >)(t|e-^ 5 >) = 5 jk 8 tr , (19) 
n o 

which is the the spanning set definition iii for the set of operators \D) = D, 
with dual (£)t|. = Tr[Z)t.]. 

Then, it is possible to write the spin tomography identity as 



A = 



2s + 1 
4tt 2 



2tt 



dn J dif) sin 2 | Tr [AD^^j, n)] D(i/>, n) , 



(20) 



a, o 

from which the following reconstruction procedure is derived 



(i) = 2 J±± Idn p(m, n) / sin 2 t Tr \A e M s ^- m ) 
47r m =~ s n o 2 



(21) 



where p(m, n) is the probability of obtaining m as the measurement result of 
S • n. This equation allows the reconstruction of arbitrary spin s expectation 
values (A), from spin measurements in all directions n. 

Numerical simulations show that the two preceding quorums are (for spin 
s — \) equivalent, namely the same number of experimental measurement data 
yield the same results and the same statistical error bars, apart from statistical 
fluctuations. In Fig. 1 a Monte Carlo comparison of the two spin reconstruction 
strategies based on the two different quorums is given. Both reconstructions 
are applied to a coherent spin state, defined as \a) = exp(aS' + — a*S-) \ — s), 
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Exp. data Exp. data 



Fig. 1. Monte Carlo comparison between continuous and discrete tomography for a 
spin s = \ system. Continuous tomography uses D(ip,n) as quorum, while discrete 
tomography uses the quorum Q = {a x ,a y ,a z ,i}. Left: Convergence of the mean 
value of (s z ) for a coherent a = 2 spin state for increasing number of experimental 
data (the theoretical value is given by the horizontal line). The circles o refer to 
continuous, the stars ★ to discrete tomography. Right: Plot of the statistical error 
bars of the graphs on the left vs experimental data. The error bars are obtained 
by dividing the experimental data into 20 statistical blocks. Notice that the two 
tomographic procedures are essentially equivalent. 

where S + ,S- are the spin lowering and raising operators and | — s) is the 
eigenvector of S z relative to the minimum eigenvalue. 

Weigert has shown [3] that another spin s quorum can be obtained by taking 
N s = (2s + l) 2 arbitrary^) directions and measuring the observables 
Qk — \nk)(nk\, which are the projectors for the eigenspace relative to the 
maximum eigenvalue s of the observables S -n k . We define a dual {Qk\ for the 
| Qk) by requiring 

(Qk\Qk') = hk' , (22) 

i.e. Eq. (11) of [3], which is just the dual set definition (6). Condition (22) 
together with the completeness of the chosen quorum, guarantee that \Qk) 
(with dual (Qk\) is a spanning set for £(H), thus allowing the tomographic 
identity 



Actually the choice of the directions is not completely arbitrary, but "almost 
[3] any choice yields a complete set of operators in C(7i). 
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N 3 

\A) = J2\Qk)(Q k \A), (23) 
k=i 

i. e. (using the notation of 



I 

k=l 



A = Y,Tr[AQ k ]Q k , (24) 



where Q h is the dual operator of Q k . The explicit form of the dual set Q h can be 
derived by a matrix inversion starting from Eq. (22) or by the Gram-Schmidt 
based procedure method given on page 4. The reconstruction procedure is, in 
this case, 



(A) = sJ2p(s,n k )Tr[AQ k 
k=i 



(25) 



where p(s, n k ) is the probability of obtaining the maximum eigenvalue s, when 
measuring S ■ n k . This allows the reconstruction of arbitrary spin operators A 
from measurements of the spin along N s fixed directions. 



4 Conclusions 



Recent Group Tomography [2] gives a general framework that allows to derive 
all the state reconstruction procedures that employ quorums which exhibit a 
group symmetry. Here we extended these results to generic state reconstruction 
procedures. In fact, we have seen how it is possible to give a characterization 
of tomographic formulas in terms of linear vector algebra on the vectors of the 
Liouville space of the system. 

A constructive method to derive new tomographic formulas has been pro- 
posed starting from the Gram-Schmidt orthogonalization procedure. At least 
in principle, it allows to calculate the quorum dual for the quantum systems 
that allow a discrete quorum. We have given some examples of the method in 
the spin domain, by re-obtaining all the known spin tomographies using linear 
vector algebra arguments. For the sake of illustrating the method, we limited 
our analysis to the description of spin systems, but all known tomographies 
can be analyzed in this framework [4]. Moreover, one may expect to employ 
the presented procedures to uncover new tomographies for quantum systems 
for which state reconstruction procedures are not presently known. 
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